REPORT  DOCUMENTATION  PAGE 


AFOSR-TR-97 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searchii 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesti 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  222D2-430l  and  to  the  Office  of  Management  and  Budget,  Paperwo 


1.  AGENCY  USE  ONLY 2.  REPORTDATE  3.  REPORT  i  iru 

August  1996  Final  Technical  Report  15  Aug  93  to  14  Aug  96 


I  3.  REPORT  I 


■<3  -.3 


4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

Characterization  of  Joint  Nonlinear  Stiffness  and  Damping  Behavior  for  Inverse  F49620-93- 1-0462 
Dynamics  of  Flexible  Articulated  Structures 


6.  AUTHOR(S) 

Assoc.  Prof.  Brad  Paden 
Thomas  A.  Trautt 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Dept  of  Mechanical  &  Environmental  Engineering 
University  of  California,  Santa  Barbara 
Santa  Barbara,  CA  93106-5070 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITQRING  AGENCY  NAME(S}  AND  ADDRESS(ES) 

AFOSR/NA 

110  Duncan  Avenue,  Rm  B115 
Bolling  AFB,  DC  20332-8050 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


F49620-93- 1-0462 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

An  inverse  dynamics  algorithm  is  derived  for  active  vibration  quenching  of  structures.  The  algorithm  uses  frequency 
domain  techniques  to  compute  an  input  function  needed  to  produce  a  desired  response  at  a  particular  degree  of  freedom.  Th 
desired  response  is  a  transition  from  the  initial  vibrating  condition  to  a  non-vibration  condition.  The  algorithm  can  also  be 
used  to  modify  the  input  function  to  correct  for  system  disturbances  while  the  input  function  is  already  being  applied  to  the 
system. 

The  algorithm  is  demostrated  in  a  simulation  of  a  simply  supported  beam  controlled  by  a  torque  actuator  at  one  end  of  th 
beam.  The  finite  element  method  is  used  to  discretizedA^fe^uations  of  motion  of  the  beam. 


^9971021  260 


14.  SUBJECT  TERMS 

^R  OF  PAGES 

10 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

UL 

Standard  Form  298  Rev.  2-89  EG 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


Final  Technical  Report 
Air  Force  Office  of  Scientific  Research 
Grant  Number  F49620-93-1-10462 


Research  Topic:  Characterization  of  Joint  Nonlinear  Stiffness  and  Damping  Behavior  for 
Inverse  Dynamics  of  Flexible  Articulated  Structures. 

Principal  Investigator:  Brad  Paden 

Student:  Thomas  A.  Trautt 


Student's  Progress  for  the  Period  August  1993  to  August  1996: 

The  student  has  completed  a  Ph.D.  in  Mechanical  Engineering  and  a  Master  of  Arts  in  Applied 
Mathematics.  The  title  of  the  dissertation  for  the  Ph.D.  is  "Inverse  Dynamics  of  Flexible 
Articulated  Structures  with  Non-Zero  Initial  Conditions."  Technical  papers  completed  and  in 
progress  are  listed  below. 


Paper  presented  at  conference: 

Trautt,  T.  A.  and  Bayo,  E.,  "Active  Vibration  Quenching  using  Inverse  Dynamics",  Proceedings 
of  the  1995  Design  Engineering  Technical  Conferences,  ASME,  DE-Vol.  84-3,  Vol.  3,  Part  C, 
PP.  239-246,  September,  1995. 


Paper  published: 

Trautt,  T.  A.  and  Bayo,  E.,  "Inverse  Dynamics  of  Non-Minimum  Phase  Systems  with  Non-Zero 
Initial  Conditions,"  Journal  of  Dynamics  and  Control,  vol.  7,  pp.  49-71,  1997. 


Paper  submitted  for  publication: 

Trautt,  T.  A.  and  Bayo,  E.,  "Active  Vibration  Quenching  of  Structures  using  Inverse 
Dynamics,"  submitted  to  the  Journal  of  Sound  and  Vibration. 


Paper  in  progress: 

Trautt,  T.  A.  and  Bayo,  E.,  "Inverse  Dynamics  of  Flexible  Manipulators  with  Coulomb  Friction 
or  Backlash  and  Non-Zero  Initial  Conditions." 


DE-Vol.  84-3, 1995  Design  Engineering  T^hnical  Conferences 

Voiume  3  —  Part  C 
ASME  1995 


ACTIVE  VIBRATION  QUENCHING  USING  INVERSE  DYNAMICS 


Thomas  A.  Trautt 

Department  of  Mechanical  Engineering 
University  of  California 
Santa  Barbara,  California 


Eduardo  Bayo 

Department  of  Mechanical  Engineering 
University  of  California 
Santa  Barbara,  California 


ABSTRACT 

An  inverse  dynamics  algorithm  is  derived  for  active 
vibration  quenching  of  structures.  The  algorithm  uses 
frequency  domain  technicques  to  compute  an  input  function 
needed  to  produce  a  desired  response  at  a  particular  degree  of 
freedom.  The  desired  response  is  a  transition  from  the  initial 
vibrating  condition  to  a  non-vibrating  condition.  The 
algorithm  can  also  be  used  to  modify  the  input  function  to 
correct  for  system  disturbances  while  the  input  function  is 
already  being  applied  to  the  system. 

The  algorithm  is  demostrated  in  a  simulation  of  a  simply 
supported  beam  controlled  by  a  torque  actuator  at  one  end  of 
the  beam.  The  finite  element  method  is  used  to  discretize  the 
equations  of  motion  of  the  beam. 


1.  INTRODUCTION 

Active  vibration  quenching  may  be  useful  for  structures 
where  passive  damping  is  not  practical  such  as  space  based 
structures.  Inverse  dynamics  may  be  applicable  to  active 
vibration  quenching  of  such  structures.  Inverse  dynamics  of 
structures  is  the  process  of  determining  load  profiles  to 
produce  a  given  displacement  profile  as  a  function  of  time  at 
some  point  of  a  structure.  Forward  dynamics  is  the  process  of 
finding  displacements  given  the  loads.  Forward  dynamics  is 
much  simpler  than  inverse  dynamics  because  it  is  causal 
whereas  inverse  dynamics  may  be  non-causal.  Because 
forward  dynamics  problems  are  causal,  they  can  be  solved  by 
integrating  forward  in  time.  Many  numerical  integration 
schemes  are  available  for  the  forward  dynamics  problem.  In 
the  non-causal  problem  encountered  in  inverse  dynamics,  the 
load  profile  to  be  computed  begins  before  the  known 
displacement  profile  begins  so  the  solution  cannot  be 
determined  by  integrating  forward  in  time.  Inverse  dynamics 


requires  methods  which  look  both  forward  and  backward  in 
time.  One  method  which  has  been  used  for  integrating  in  both 
directions  in  time  is  to  use  frequency  domain  techniques 
involving  Fourier  transforms.  Another  method  is  to  use 
convolution  integrals  in  the  time  domain.  A  third  method  is  to 
separate  the  system  into  causal  and  non-causal  components 
and  integrate  the  non-causal  component  backwards  in  time  and 
the  causal  component  forward  in  time. 

Most  of  the  research  in  inverse  dynamics  has  been  for 
applications  in  controling  flexible  robot  arms.  Bayo  (1987) 
introduced  the  inverse  dynamics  method  with  frequency  domain 
techniques  for  open  loop  control  of  a  single-link  flexible 
robot.  The  position  of  the  tip  of  the  flexible  link  was 
controlled  by  a  torque  actuator  at  the  pinned  end  of  the  link. 
The  problem  was  solved  for  the  link  initially  at  rest.  The 
method  was  validated  experimentally  by  Bayo,  Movaghar,  and 
Medus  (1988).  A  similar  method  presented  by  Bayo  and 
Moulin  (1989)  was  solved  in  the  time  domain  by  converting 
the  Fourier  transform  equation  into  a  convolution  integral. 

Chapnik,  Heppler,  and  Aplevich  (1993)  used  inverse 
dynamics  with  frequency  domain  techniques  to  control  the 
impact  response  of  a  flexible  link. 

Lopez-linares  et  al.  (1991)  used  inverse  dynamics  with 
frequency  domain  techniques  and  feedback  to  control  a  flexible 
link.  Inverse  dynamics  was  used  to  compute  feed  forward 
torques  with  the  link  initially  at  rest  and  feedback  control  was 
used  to  correct  for  initial  conditions  and  disturbances.  Kwon 
and  Book  (1990)  presented  a  similar  approach  using  a 
combination  of  inverse  dynamics  and  feedback  control.  The 
inverse  dynamics  was  computed  by  separating  the  system 
equations  into  causal  and  non-causal  components. 

Inverse  dynamics  has  also  been  used  for  multibody  flexible 
manipulators.  Bayo,  Papadopoulos,  and  Stubbe  (1989)  and 
Ledesma  and  Bayo  (1993)  solved  multibody  problems  with  the 
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system  initially  at  rest.  Paden  et  al.  (1993)  used  inverse 
dynamics  and  feedback  to  control  a  multibody  flexible 
manipulator.  The  inverse  dynamics  was  computed  for  the 
system  initially  at  rest  and  the  feedback  control  was  used  to 
correct  the  trajectories  for  initial  conditions  and  disturbances. 

Jayasuriya  and  Choura  (1990)  presented  a  method  for 
vibration  quenching  which  was  demonstrated  on  a  pinned- 
pinned  beam  with  a  single  fixed  point  actuator.  The  forcing 
function  was  a  sum  of  sines  and  cosines  at  the  natural 
frequencies  of  the  system.  The  coefficients  to  the  sines  and 
cosines  were  determined  so  that  all  the  displacements  and 
velocities  were  zero  at  some  future  time.  This  method  controls 
the  structure  at  the  end  point  of  a  time  interval  as  in  the 
slewing  problem  whereas  the  inverse  dynamics  methods 
mentioned  above  control  a  degree  of  freedom  through  a 
predetermined  trajectory  during  a  time  interval. 

In  the  present  paper,  inverse  dynamics  is  used  to  quench 
vibrations  in  a  pinned-pinned  beam.  A  torque  actuator  acts  on 
one  end  of  the  beam  to  control  the  angular  deflection  at  the 
other  end  of  the  beam.  The  beam  is  modeled  with  the  same 
physical  properties  as  modeled  by  Jayasuriya  and  Choura.  The 
equations  of  motion  are  discretized  using  finite  element 
methods  assuming  Bemouli-Euler  beam  theory.  An  inverse 
dynamics  algorithm  is  derived  using  frequency  domain 
techniques.  The  algorithm  is  demonstrated  in  computer 
simulations  of  vibration  quenching  of  the  pinned-pinned 
beam.  Three  problems  are  demonstrated.  The  first  problem  is 
to  quench  the  beam  when  only  the  first  vibration  mode  is 
active.  The  second  problem  is  to  quench  the  beam  when  many 
modes  are  active.  In  the  third  problem,  a  disturbance  is  added 
to  the  system  after  the  forcing  function  is  initiated,  and  the 
forcing  function  is  then  changed  to  compensate  for  the 
disturbance.  The  ability  of  the  algorithm  to  correct  the  forcing 
function  may  lead  to  the  ability  for  real  time  control  by 
continually  updating  the  forcing  function. 


2.  MODELING 

The  inverse  dynamics  algorithm  will  be  demonstrated  in  a 
simulation  of  a  pinned-pinned  beam  as  shown  in  Figure  1.  A 
torque,  fit),  as  a  function  of  time,  t,  acts  on  the  left  end  of  the 
beam  to  control  the  angular  displacement  at  the  right  end  of 
the  beam.  The  beam  has  a  total  length  L,  Young’s  modulus  E, 
constant  cross  section  area  A,  constant  area  moment  of  inertia 
I,  and  density  p. 

We  will  use  finite  element  methods  to  discretize  the 
equations  of  motion.  We  assume  Bemoulli-Euler  beam  theory 
and  small  deflections,  and  neglect  gravity,  rotory  inertia,  and 
axial  loading.  A  typical  element  is  shown  in  Figure  2.  Let 
v(sj)  be  the  tranverse  displacement  as  a  function  of  time,  /, 
and  distance,  s,  along  the  element.  Let  v  be  the  nodal 
displacements  and  u  be  the  corresponding  nodal  loads  applied 
to  the  element.  The  potential  energy  of  the  element  is 


where  the  primes  denote  partial  derivatives  with  respect  to  s 
and  Le  is  the  length  of  the  element.  The  virtual  work  done  by 
inertia  forces  is 

fU 


inertia^  f  5v[-p  A v( 

Jo 


(2) 

where  the  dots  denote  partial  derivatives  with  respect  to  time. 
The  virtual  woric  done  by  nonconservative  forces  is 

^nc  =  (3) 


The  principle  of  virtual  work  can  be  stated 

inertia  +  —  SW  nc  —  0 


(4) 


Substituting  equations  (1),  (2),  and  (3)  into  this  equation,  we 
obtain 


pAdvv(s,t)ds  + 

.0  Jo 


pA^vv’lf  ,/)ds  +  El Sv''v''(s,t)ds  -  dv^a  =  0 


We  discretize  the  displacements  as  follows. 
v(j,f)  =  h(s)v(f)  =  v'^'h^ 


(5) 

(6) 


(7) 


In  this  equation,  h  is  a  row  vector  of  shape  functions  based  on 
the  Hermite  polynomials  (Meirovitch,  1990).  The  shape 
functions  are 

hi(s)  =  L,  fyU  -  2(jt/LeP  +  (x/L,)^] 

A2W  =  Le[ -(x/ur +  (;«/L,r] 

/•3W=l-3(x/Le)^  +  2(x/Lef 
h4(s)  =  3(xa^f-2ixR^f 

From  equation  (6)  we  can  derive  the  following  expressions, 
hv  ^ 

5v'=  5v  (h") 
v"  =  h'v 

Substituting  equations  (8)  into  equation  (5),  we  obtain 


(8) 


pA  iJhdsy  +  e  (h")Vjjv  -  u 

.  Jo  Jo 


=  0 


Setting  the  part  in  the  brackets  equal  to  zero,  we  obtain 

MeV  +  KeV  =  U 

where 

rU 

M,  =  pA  h'^hdi 


(9) 

(10) 

(11) 

(12) 

The  matrices  M,  and  K,  are  the  mass  and  stiffness  matrices 
respectively  for  an  element.  The  beam  is  modeled  with  several 
elements  as  shown  in  Figure  3.  Assembling  the  elements  we 
obtain  the  matrix  equation  of  motion  for  the  system, 


and 


Jo 

Le 


/•Le 

K«=EI  (h") 
.0 


bds 
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My  +  Ky  =  bf(/) 


(13) 


where  M  is  the  system  mass  matrix,  K  is  the  system  stiffness 
matrix,  y  is  a  vector  of  degrees  of  freedom,  and  b  is  a  loading 
vector  with  a  unit  value  at  the  degree  of  freedom  at  which  the 
torque  is  applied  and  zeros  in  the  rest  of  the  vector. 

We  add  Rayleigh  damping  to  the  system  using  the 
expression 

C=aM+PK 

where  C  is  the  damping  matrix,  and  a  and  P  are  constants.  The 
matrix  equation  of  motion  of  the  system  with  damping 
becomes 

My  +  Cy  +  Ky  =  b/(0  (15) 

The  damping  ratio  at  each  natural  frequency  can  be  calculated 
from 

?i  =  +  PtOi) 

2\coi  I  (16) 

where  is  the  damping  ratio  and  (Oj  is  the  natural  frequency  in 
radians  per  unit  time. 


3.  INVERSE  DYNAMICS 

We  want  to  solve  for  the  torque  profile,  /(/),  which  will 
quench  the  initial  vibrations  in  the  system.  Taking  the  Fourier 
transform  of  equation  (15),  we  obtain 

+  iQ>C  +  K)y(tit))=  b^co)  (17) 

or 

Ay  =b/ 
where 


A(£u)  =  (.Q)^M+i(yC  +  K) 


(18) 

(19) 


The  caret  denotes  the  Fourier  transform  and  (o  is  frequency  in 
radians  per  unit  time.  Solving  equation  (18)  for  we  obtain 
y  =  g/  (20) 


where 


g  =  A  b  +icuC  +  K)'‘b 


(21) 


Without  damping  in  the  system,  matrix  A  would  not  be 
invertible  when  (O  is  equal  to  a  natural  frequency.  For  this 
reason,  it  is  necessary  to  have  some  damping  in  the  system  for 
the  algorithm  to  be  usable. 

We  can  write  equation  (20)  in  an  expanded  form  as  follows. 


'y  1^ 

iy<i  1 

\gn 

/ 


(22) 


The  kth  equation  corresponding  to  the  controlled  degree  of 
freedom,  gives  us 


yk  =  gkf 

A 

Solving  for  /  ,  we  obtain 
/  =yk  i  gk 


(23) 


(24) 


We  need  to  define  a  profile  for  y^(/)  which  will  cancel  the 
vibrations  due  to  the  initial  conditions.  We  also  need  the 
profile  to  be  finite  in  time  so  that  we  can  compute  its  fast 
Fourier  transform  (FFT).  To  cancel  the  vibrations  for  a  finite 
time,  we  need  the  profile  to  be  the  negative  of  the 
homogeneous  response  for  a  finite  time.  Therefore,  we  create  a 
windowed  profile,  z(t),  as  follows. 


z{t)  -  -w(r)  yk{t) 


(25) 


where  w(t)  is  a  windowing  function  and  y^(t)  is  the 
homogeneous  response  to  the  initial  conditions.  The 
homogeneous  response  can  be  computed  using  numerical 
integration.  The  general  shape  of  the  windowing  function  and 
its  first  two  derivatives  are  shown  in  Figure  4.  The  time  from 
to  to  tf  is  the  window  of  time  for  the  FFT  computation.  We  will 
go  into  more  detail  on  the  windowing  function  later  in  our 

discussion.  We  compute  the  FFT  of  z(t)  and  substitute  it  for  y* 
in  equation  (24)  to  obtain 


f-z/gk 


(26) 


where  gk  is  computed  at  the  same  frequencies  used  to  compute 
the  FFT.  Taking  the  inverse  fast  Fourier  transform  (IFFT)  of 

A 

f(co)y  we  obtain  the  torque  profile,  yfr).  We  can  think  of  fit)  as 
a  composition  of  two  parts.  The  first  part  cancels  the 
vibration  and  the  second  part  causes  the  system  to  transition 
back  to  the  vibrating  condition.  We  need  to  eliminate  the  part 
of  the  forcing  function  which  causes  the  system  to  transition 
back  to  the  vibrating  condition.  We  do  this  by  eliminating 
the  second  half  of  fit)  as  follows. 


^2(0  = 


\fit)  ,  to^^tir 

1  0  »  tm  <  ^ 


(27) 


In  this  formula,  to  is  the  starting  time  for  the  FFT  window, 
is  the  midpoint  of  the  FFT  window,  and  f2(t)  is  the  desired 
forcing  function  which  will  cancel  the  vibrations.  There  will 
be  some  error  due  to  truncation  and  overlap  of  the  two  parts  of 
the  torque  profile,  fit).  This  error  can  be  made  negligible  by 
choosing  large  enough  time  spans  in  the  windowing  function. 
In  the  time  spans  where  z(t)  is  following  a  homogeneous 
response,  the  torque  profile  approaches  zero  as  the  rest  of  the 
system  approaches  the  homogeneous  response. 

We  may  want  to  modify  the  forcing  function  while  it  is 
being  applied  to  correct  for  disturbances  in  the  system.  We 
assume  we  already  have  a  forcing  function,  fj(t),  which  will 
produce  a  profile  of  which  will  follow  a  homogenous 
response  after  a  time  I2  and  has  zero  force  after  t^.  The  forcing 
function. /2(r)  from  equation  (27)  will  have  these  properties. 
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Using  numerical  integration,  we  compute  the  response  of  y^t) 
due  to  initial  conditions  and  the  forcing  function  from 

My  +  Cy  +  Ky  =  b/i(r)  ^28) 

The  windowed  profile  is  then  computed  from 

z(r)  =  -w(r)yi(r)  ^29) 

We  compute  giica)  at  discrete  frequencies  using  the  following 
equation. 

g  =  (-tu^M  +itoC  +  K)‘'b  (3QJ 

We  compute  the  FFT  of  z(t)  and  then  solve  for  /  from 

/  =  W9*  ^  (31) 

We  compute  the  IFFT  of  /(<o)  to  obtain  f(t)  and  add  the  first 
half  of  f(t)  tofjft)  to  obtain  the  new  desired  forcing  function 


f2(.t)  = 


,  to  it  Stm 

0  ,  tm<  / 


Using  numerical  integration  we  can  compute  the  new  response 
of  the  system  from 

My  +  Cy+  Ky  =  b/2(r)  (33^ 

Equations  (28)  through  (33)  summarize  the  algorithm.  The 
algorithm  can  be  continually  repeated  by  using  the  profile  of 
f2(0  or  a  portion  of/2(r)  with  an  apropriate  time  shift 
substituted  for  fj(t)  in  equation  (28). 

The  derivation  of  the  windowing  function,  w(t),  follows. 
We  want  w(t)  to  produce  smooth  functions  of  z,  z  and  z  so 
that  the  resulting  torque  profile,  fit),  is  also  smooth.  For  z  to 
be  smooth,  we  need  z  *  to  be  continuous.  The  windowed  profile 
and  its  derivatives  are 

z--wyk 

Z  =  -w  yk  -w  yk 

z  =  -ivyk  •  2w  yk  ‘  wyk 

Z  =  -wyk  -  3wyk  -  3Wyi  -  w  yi  ^3^^ 

The  windowing  function,  w(/).  transitions  from  zero  to  one 
between  tj  and  t2.  For  z,  z  and  z*  to  be  smooth  at  t]  and  t2,  we 
have  the  following  requirements  on  w(t). 

Mti)  =w(ti)  =w(ti)  =w'(ti)  =  0 
w(t2)=l 

H'(t2)  =vv(t2)  =vv*(t2)  =  0 

T.  .  (35) 

Between  times  tj  and  t2,  we  let  w(t)  be  described  by  a 

polynomial,  p(9) . 

^(t)=p{e) 

where 

t2  -  ti  (37) 


Then,  p( 6)  will  have  the  following  requirements. 

P(0)  =  p'(0)  =  p"(0)  =  p"'(0)  =0 
/>(!)  =  1 

P'(1)=P"(1)  =/>"'(!)  =  0  ^3g^ 

The  primes  denote  derivatives  with  respect  to  B.  We  have 
eight  constraints  on  p(0)  so  we  need  a  polynomial  with  eight 

constants.  The  polynomial,  p(d)  ,  and  its  derivatives  can  be 
written  as 

p(8)  =  ao+aie+a20^+a30^+a40\a50^+a6e*+a70’ 
p'(6)  =  ai+2a20+3a30^+4a40^+5a50V6a60^+7a70* 
p'\9)  =  2a2+6a3^12a40^+2Oa5fl^+3Oa6flV42a70^ 
p"\9)  =  6a3+24a40+6Oa50^+12Oa60^+21Oa7e‘‘  (39) 

where  ao  through  are  constants.  From  the  requirements 
p(0)=pX0)=/>"(0)=p'"(0)=0,  we  obtain 

ao=a,=a2  =  a3  =  0  ^40^ 

From  the  other  four  requirements  we  obtain 

34  +  a5  +  a6  +  a7  =  l 

4a4  +  Sas  +  6a6  +  7a7  =  0 

12a4  +  20a5  +  30a6  +  42a7  =  0 

24a4  +  eOas  +  12036  +  210a7  =  0  f4n 


or  in  matrix  form 

'1111  iraal  n 
4  5  6  7  aj  _  0 

12  20  30  42  36  0 

-  24  60  120  210  JLa7J  L  0 


The  solution  to  these  equations  is 


a4']  35 

as  ^  -84 

ae  70 

a7J  L  -20 


and  we  obtain  the  polynomial 


pfS)  =  -200’  +  700*  -  840*  +  350“* 


where 

t2-ti 


tl  ^/^t2 


Between  13  and  t4,  we  can  use  p( ^)  where  (p  is  given  by  the 
following  expression  which  will  give  us  a  transition  from  one 
back  to  zero. 

U-t3  .  t3^/^t4  (46) 
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Putting  the  pieces  of  w(t) 
formula. 

together,  we 

obtain  the  following 

'  0 , 

VI 

VI 

S 

I 

Pirn) . 

tl  ^  f  <  t2 

w(0=  i 

1 , 

t2<t  <t3 

p(m . 

t3<f  <4 

>  0, 

t4<t  <tf 

(47) 

This  formula  was  used  to  make  the  plots  of  w,  vv,  and  vv  shown 
in  Figure  4. 


4.  APPLICATION 

The  inverse  dynamics  algorithm  will  be  demonstrated  in  a 
computer  simulation  of  a  pinned-pinned  beam  as  shown  in 
Figure  1.  The  beam  is  modeled  with  six  finite  elements  and  has 
the  following  properties: 

L  =  0.762  m  =  total  length 

A  ss  1.613  X  10"^  m^  =  cross  sectional  area 

I  =  9.6361  X  10**^  m"^  =  area  moment  of  inertia 

p  =  2710  kg/m^  =  density 

E  =  71  X  10^  N/m^  =  Young’s  modulus 

a  =  .01/sec,  p  =  .0001  sec,  Rayleigh  damping  constants 

(48) 

The  windowing  function,  wfr),  shown  in  Figure  4  will  have 
the  following  values  for  to  through  tf. 

to  =  0.0  sec 
t|  =0,1  sec 
t2  =  0.2  sec 
tjn  =  0.3  sec 
t3  =  0.4  sec 
t4  =  0.5  sec 

tf  =  0.6  sec  (49) 

The  first  six  natural  frequencies  of  the  discretized  system  are 

0)1  =  3.38  Hz  =  21.2  rad/scc 
0)2  =  13.5  Hz  =  84.8  rad/sec 
0)3  =  30.6  Hz  =  192  rad/sec 
0)4  =  54.7  Hz  =  344  rad/sec 
0)5  =  86.8  Hz  =  545  rad/sec 

0)5  =  135  Hz  =  848  rad/sec  (50) 

The  damping  ratios  associated  with  the  first  six  natural 
frequencies  are 

=  .0013 
^2  =  .0043 


^3  =  .0096 
C4  =  .017 
Cs  =  .027 

':6  =  .042  (51) 

We  can  see  that  this  system  is  very  lightly  damped  at  the  lower 
frequencies. 

The  inverse  dynamics  algorithm  will  be  demonstrated  on 
three  problems.  The  first  problem  is  to  quench  the  vibrations 
when  only  the  first  mode  is  active.  The  second  problem  is  to 
quench  the  vibrations  when  many  modes  are  active.  In  the 
third  problem,  a  disturbance  is  added  to  the  system  while  the 
torque  profile  is  already  being  applied,  and  the  inverse 
dynamics  algorithm  is  used  to  change  the  torque  profile  to 
correct  for  the  disturbance. 


ill — First  Mode  Active 

The  initial  displacement  profile  of  the  beam  is  the  first  mode 
shape  with  a  displacement  of  .075  m  at  the  center  of  the  beam. 
The  initial  velocities  are  all  zero.  The  inverse  dynamics 
algorithm  is  used  to  compute  a  torque  profile  which  will  cancel 
the  vibrations  in  the  beam.  The  response  of  the  beam  is 
computed  using  Newmark's  constant-average-acceleration 
method  (Bathe,  1982).  The  torque  profile  and  displacement 
profiles  are  shown  in  Figure  5.  Angular  displacements  at  the 
right  end  of  the  beam  and  mid-span  displacements  are  shown 
with  and  without  the  torque. 

We  see  that  the  torque  is  at  a  negligible  level  at  0.0  second 
and  at  0.3  second  where  the  torque  profile  is  truncated.  These 
times  correspond  with  to  and  shown  in  the  graph  of  the 
windowing  function  in  Figure  4. 

The  angular  dispacement  which  is  the  controlled  degree  of 
freedom  settles  down  by  0.2  second  which  corresponds  with 
the  time  t2  shown  in  the  graph  of  the  windowing  function. 
The  torque  profile  continues  to  damp  the  rest  of  the  system  for 
a  short  time  after  the  controlled  degree  of  freedom  has  already 
settled  down.  We  see  that  the  mid-span  displacement  settles 
down  after  the  angular  displacement  settles  down. 

— Many  Modes  Active 

The  initial  conditions  of  the  vectors  y  and  y  are  all  zero  except 
for  an  angular  velocity  of  10  rad/sec  at  the  left  end  of  the  beam. 
These  initial  conditions  include  many  vibration  modes.  The 
inverse  dynamics  algorithm  is  used  to  compute  the  torque  to 
quench  the  vibrations.  The  torque  profile  and  displacement 
profiles  with  and  without  the  torque  arc  shown  in  Figure  6.  As 
in  the  previous  problem,  we  see  that  the  angular  displacement 
settles  down  by  0.2  second  and  the  mid-span  displacement 
settles  down  after  the  angular  displacement  settles  down. 

iJi — Correct  for  a  Disturbance 

Starting  with  the  previous  solution,  we  add  a  disturbance  to 
the  torque  profile  from  time  t=0.02  second  to  time  t=0.03 
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second.  At  time  t=0.06  second,  we  apply  the  inverse  dynamics 
algorithm  to  correct  the  torque  profile.  Torque  profiles  with 
the  disturbance  and  displacement  profiles  are  shown  in  Figure 
7  with  and  without  the  correction. 

Since  the  algorithm  to  correct  the  profile  is  started  at  time 
t=0.06,  we  expect  the  angular  displacement  to  settle  down  by 
0.26  second.  As  we  see  in  the  plot,  the  angular  displacement 
does  settle  down  by  0.26  second  as  expected. 


5.0  CONCLUSION 

An  inverse  dynamics  algorithm  for  quenching  vibrations 
has  been  presented.  The  algorithm  was  demonstrated 
analytically  on  a  pinned-pinned  beam  with  a  torque  actuator  at 
one  end  and  has  shown  good  results.  The  vibrations  are 
damped  out  quickly  and  the  required  torque  profiles  are  smooth. 
The  algorithm  can  be  repeated  to  correct  for  errors  and 
disturbances  after  the  forcing  function  is  initiated.  This  will 
lead  to  the  ability  for  real  time  control  by  continually  updating 
the  algorithm  at  discrete  time  intervals. 

The  algorithm  can  be  applied  to  more  complicated 
structures.  The  torque  or  force  actuator  should  be  applied  to  a 
degree  of  freedom  which  is  not  a  node  of  one  of  the  natural 
frequencies.  The  actuator  must  also  be  non-colocated  with  the 
controlled  degree  of  freedom.  The  controlled  degree  of  freedom 
should  not  be  at  a  node  of  any  natural  frequency  which  is  to  be 
quenched. 

The  algorithm  is  robust  to  errors  in  initial  conditions. 
Because  the  algorithm  uses  superposition,  we  can  superpose 
the  assumed  initial  conditions  and  the  error  in  the  initial 
conditions.  Applying  the  algorithm  will  eliminate  the 
vibrations  due  to  the  assumed  initial  conditions.  The  response 
after  the  algorithm  is  applied  will  be  the  homogeneous 
response  to  the  error  in  the  initial  conditions.  Each  time  the 
algorithm  is  applied  this  error  will  be  reduced. 

The  algorithm  is  also  robust  to  errors  in  damping,  but  is  not 
very  robust  to  errors  in  mass  and  stiffness.  Errors  in  mass  and 
stiffness  can  cause  the  computed  response  to  be  out  of  phase 
with  the  actual  response  especially  at  the  higher  frequencies 
and  we  would  not  get  the  desired  cancelling  effect  when 
applying  the  algorithm.  The  robustness  can  be  improved  by 
using  inverse  dynamics  to  compute  a  feed  forward  control  and 
adding  closed-loop  tracking  control  to  track  the  trajectories 
computed  by  the  inverse  dynamics  algorithm  (Kwon  and  Book, 
1990),  (Lopez-Linares  et  al.,  1991  )  and  (Paden  et  al.,  1993). 

More  research  is  needed  to  make  the  computations  more 
practical  for  real  applications.  Simpler  models  could  be  used  if 
only  the  first  one  or  two  vibration  modes  need  to  be  actively 
quenched.  If  we  assume  only  the  first  two  modes  are  present, 
then  less  sensors  will  be  needed  to  estimate  the  initial 
conditions  of  the  system.  Simplifying  the  model  and  reducing 
the  number  of  sensors  would  reduce  the  computations  at  the 
cost  of  losing  some  accuracy.  The  computations  can  also  be 
made  more  efficient  by  using  the  time  domain  approach 
presented  by  Bayo  and  Moulin  (1989)  instead  of  the  frequency 
domain  approach. 
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FIG.  1.  PINNED-PINNED  BEAM 
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FIG.  2.  BEAM  ELEMENT 


FIG.  3.  FINITE  ELEMENT  MODEL 


Time,  (sec) 


FIG.  4.  WINDOWING  FUNCTION,  w(t)  (SOLID  LINE), 
AND  ITS  DERIVATIVES,  vi-CO  (DASHED  LINE)  AND 
w(f)  (DOTTED  LINE),  NORMALIZED  TO  ONE. 


FIG.  5.  RRST  MODE  ACTIVE  -  TORQUE  TO  CANCEL 
VIBRATIONS,  ANGULAR  DISPLACEMENT  AT  RIGHT  END, 
MID-SPAN  DISPLACEMENT,  WITH  TORQUE  (SOLID 
LINE),  WITHOUT  TORQUE  (DASHED  LINE). 


V 


FIG.  6.  MANY  MODES  ACTIVE  -  TORQUE  TO  CANCEL 

VIBRATIONS,  ANGULAR  DISPLACEMENT  AT  RIGHT  END  DISTURBANCE  IN  TORQUE  PROFILE  -  TORQUE . 

MID-SPAN  DISPLACEMENT.  WITH  TORQUE  (SOLID  ’  ANGULAR  DISPLACEMENT  AT  RIGHT  END,  MID-SPAN 
LINE),  WITHOUT  TORQUE  (DASHED  LINB  DISPLACEMENT,  WITH  CORRECTION  (SOLID  LINE). 

''  WITHOUT  CORRECTION  (DASHED  LINE). 
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